Controlling the Sign of Magnetoconductance in Andreev Quantum Dots 
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We construct a theory of coherent transport through a baUistic quantum dot coupled to a super- 
conductor. We show that the leading-order quantum correction to the two-terminal conductance 
of these Andreev quantum dots may change sign depending on (i) the number of channels carried 
by the normal leads or (ii) the magnetic flux threading the dot. In contrast, spin-orbit interaction 
may affect the magnitude of the correction, but not always its sign. Experimental signatures of the 
effect include a non-monotonic magnetoconductance curve and a transition from an insulator-like 
to a metal-like temperature dependence of the conductance. Our results are applicable to ballistic 
or disordered dots. 

PACS numbers: 74.45.-l-c, 74.78.Na, 73.23.-b 



Introduction. Low temperature experiments on dif- 
fusive metals coupled to superconductors have reported 
large interference effects analogous to coherent backscat- 
tering, weak-localization and Aharonov-Bohm oscilla- 
tions [H, [3, [3, 0, H, H, 0], one to two orders of magni- 
tude above the universal amplitude 0{e]Jh) they have in 
purely metallic mesoscopic conductors [8[. In some cases, 
a weak localization-like behavior, in the form ofpositive 
magnetoconductance near zero field is observedl2U7|, but 
often one sees negative magnetoconductance [JU a Q . 
Theoretical works predict that Andreev reflection from 
the superconductor induces this large quantum correc- 
tion to transport 0, • The general expectation is that 
this correction is similar to a magnified weak-localization 
correction, in that its sign is determined bythe presence 
or absence of spin-orbit interaction (SOI) [8|, lll|. In this 
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paper we revisit this issue, and find that this interfer- 
ence correction has very different properties from weak- 
localization. In particular, we show that both the specific 
lead-geometry and an applied magnetic flux can reverse 
its sign, while SOI need not. 



Andreev reflection 1^ is the dominant low energy pro- 
cess at the interface between a metal and a superconduc- 
tor. It involves an electron (hole) being retroreflected 
as a hole (electron) and retracing the path previously 
followed by the electron (hole). In this article, we ex- 
tend the trajectory-based semiclassical theory to include 
Andreev reflection, analyze the conductance of a two- 
dimensional ballistic quantum dot coupled to one super- 
conducting lead and two normal leads, as in Fig. [T] We 
dub this system an Andreev quantum dot. We arrive at 
the surprising conclusion that the interference effects can 
be reversed from localizing to antilocalizing by changing 
the widths of the normal leads, or by threading a frac- 
tion of a magnetic flux quantum through the dot. In 
contrast SOI need not cause such a reversal. This is 
very different from weak-localization in purely metallic 
conductors, whose sign is solely determined by the pres- 
ence or absence of SOI [1, [U, [l^l ■ We predict two clear 
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Figure 1: (a) A two-dimensional Andreev quantum dot in a 
three-terminal geometry, with two normal (N) and one su- 
perconducting (S) lead. (b,c) The two possible two-terminal 
set-ups obtained from such a dot. Either (b) the S lead is 
contacted to one of the N leads, or (c) the S lead is floating. 



experimental signatures of these interference effects in 
the form of non-monotonic magnetoconductance curves 
(see Fig. [3]) and a transition from an insulator-like to 
a metal-like temperature-dependence of the conductance 
as one changes the magnetic field or the ratio of the lead 
widths. This transition occurs because thermal averaging 
destroys quantum interferences, thus depending on the 
sign of the effect, the conductance increases or decreases 
by many times /h as the temperature increases. 

Semiclassical transport with superconductivity. 
According to the scattering approach to transport, the 
current in normal lead i is given by [l^ 
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Figure 2: Contributions to (T^^) (first three) and (T^j^) (last four) considered in this letter. Thick green (thin violet) paths 
indicate electrons (holes), dashed lines indicate complex-conjugated amplitudes. Normal leads are labelled i,j while the S lead 
is superconducting. The contributions are classified by the number of uncorrelated Andreev reflections (eeO has none, ee2I and 
ee2II both have two). The full (open) squares on the S lead indicate a factor of rj (77*) and the ellipses mark encounters. 



where fj,j is the chemical potential of normal (N) lead j 
and fj-s of all the superconducting (S) leads. The Fermi- 
Dirac distribution, /(e), has e measured from ns- Here 
we use trajectory-based semiclassics to find the scatter- 

e, h 
fix- 



ing probability T^^ to go from quasiparticle /3 
(electron, hole) in lead j to quasiparticle a in lead 



tending trajectory-based semiclassics Il5i 
to include Andreev reflection, one has |2( 
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This expression sums over all classical trajectories 7I and 
72 entering the cavity at yo on a cross-section of lead j 
and exiting at ?/ on a cross-section of lead i, while convert- 
ing a P quasiparticle into an a quasiparticle. The phase 
SS = S^i — Sj2 gives the difference in action phase ac- 
cumulated along 7I and 72, and gives the stability of 
the trajectory 7. In contrast to Ref. [20|, we consider the 
physically more prevalent situation of an Ehrenfest time 
negligible against the dwell time td inside the dot. In 
that case, even with perfect Andreev reflection, quantum 
uncertainties combined with the chaotic dynamics make 
the retroreflected quasiparticle diverge from the incom- 
ing quasiparticle path well before it leaves the dot 21 1. 



Therefore classical paths undergoing Andreev reflections 
consist of electron and hole segments that do not neces- 
sarily retrace each other all the way. For transmission 
probabilities {T°^^) averaged over energy or dot shape, 
we must pair the paths 7I and 72 in Eq. ([2]) in ways 
that render their action phase difference stationary. To 
do this we either pair a path with a complex conjugate 
path, e-e* or h-h*, or we pair an electron path with a hole 
path, e-h or e*-h*. Path-pairs can meet and swap pair- 
ings at encounters^ as shown in Fig.[2l Following Ref. [l3| 
we distinguish between encounters that lie entirely inside 
the dot and those that touch a lead. 

Feynman rules. Contributions relevant to current 
noise in purely metallic samples H, 17, 18, 31, become 
relevant for the calculation of the current itself in the 
presence of S leads when they can be made from only 
two classical trajectories with some segme nts as electron 
and others as holes. FromRefs. [l3,[llllll and the above 
considerations, we derive the following Feynman rules for 



calculating transmission through an Andreev quantum 
dot. The dot is connected to normal and superconducting 
leads, each carrying iV^ ^ 1 and N^j ^ 1 transport chan- 
nels respectively, and we write iVx = X^i + -^s^ 
For a perpendicular magnetic field, b = B/Bc, measured 
in units of the field Be — (/i/ej4)(ro/TD)^^^ that breaks 
time-reversal (TR) symmetry in a quantum dot of area 
A with time of flight tq, the Feynman rules read: 

(i) An c-c* or h-h* path-pair gives a factor of [Nt{1 + 
x6^)]^^, with X = 1 for time-reversed paths and x = 
otherwise. 

(ii) An e-h or e*-h* path-pair gives N,^^ x (1 ± i2eTD + 
X6^)~^, with upper (lower) sign for e-h (e*-h*). 

(iii) An encounter inside the dot and connecting e,e*, h 
and h* paths (as in he2II) gives a factor ~Nt- 

(iv) An encounter inside the dot and connecting e, e, 
c* and h paths (as in ee2II) gives a factor of — A^t(1 + 
i2eTD -I- 6^); this factor is complex conjugated (c.c.) for 
an encounter connecting e, e*, e* and h* paths. 

(v) An encounter touching a N lead i (S lead j) gives a 
factor of A^i (Nsj)- 

(vi) A path-pair that ends at a N lead i (S lead j), while 
not in an encounter, gives a factor of Ni (Nsj). 

(vii) Andreev reflections at S leads involving e— >h give 
a factor of jye"'*^-! while those involving h— >e give a 
factor of ?7e'*^j (e* ^h* and h* — >e* give the c.c. of 
these factors), where <I>sj is the S phase on lead j, and 
1] ~ cxp[— iarcos (e/A)] is the Andreev reflection phase. 
We note that these rules equally follow from random- 
matrix theory [l7l |. 

In our analysis of the consequences of these rules, we 
consider temperatures well below the superconducting 
gap A where = — i, and consider a single S lead (set- 
ting (f>s = without loss of generality). The rules indi- 
cate that a path-pair going from encounter to encounter 
reduces the contribution by a factor of ©[A't]- Thus to 
leading order in iVr, we can neglect such (weak localiza- 
tion) contributions. This does not restrict the number of 
encounters, because the price to add an encounter whose 
additional legs go to S leads is 0[(7Vs/^t)^]- Wc there- 
fore take Ns/Nt <^ 1, and expand in the number of 
uncorrelated Andreev reflections. 

Restricting ourselves to 0[{Ns/Nt)'^], we need to con- 
sider the contributions shown in Fig. [2] involving no more 
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than two uncorrelated Andreev reflections. The contri- 
butions to are 



(7^7°) = N,N,/Nt, 
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The contributions to (T'/j") are 
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Semiclassics gives {Tl^) = (T^f ) and (T;f ) = (T;!j"=). 
These contributions preserve unitarity up to and includ- 

ingO[(7Vs/E,^^.)']- 

Set-up v^rith an S lead. We first consider the sit- 
uation where the S lead's potential is fixed externally. 
This may be the three-terminal device of Fig. [1^ with 
both the R and S leads grounded, while the L lead is 
biased at electrochemical potential /xl = eV . Alterna- 
tively this may be the two-terminal device in Fig. [TJa 
with the S and R leads join at a bulk contact (with con- 
tact conductance vastly greater than the dot), a macro- 
scopic distance away from the dot. In either case, the L 
lead current is /l = (2e^//i) [^ci + (^^qmCT, 6)] F, where 
we define a dimensionless classical conductance gc\ = 
iVL(iVR + 2Ns)/{Ni^ + iVR + 2Ns) For Ns <. N 

the quantum interference correction is 



iVL[jVR-4iVL(l + 6^)]jV| ^^y 



Sg^i. (5) 



In the regime of experimental interest the weak- 
localization correction in the absence of the S lead 
Sg^i ~ -AlA^r/[(A^l + A^r)^(1 + b^)] is smah enough 
to neglect (as in Fig. [3]). The e- integral in Eq. ^ 
with {T"/) in Eqs. (HI) leads to f{T,b) = aC{2,l/2 + 
(1 -I- b'^)a), with a ~ / 'iiTk-QT for a Thouless energy 
-E-Th = fi/TB, and the generalized ^-function C{2,x) = 
t exTp[—xt]/ (1 — exp[— t])dt. This gives the two asymp- 
totics /(T ~> oo,6) 7r£;T/[8fcBr(l -I- b^)] and f{T ^ 
0,b) ^ l/(l + &2)2. 

At zero temperature, we find three regimes for Sgq^i- 

(a) For Nr < 2Ni^, Sgq^n < for all values of 6, and gives 
a monotonic magnetoconductance curve. 

(b) For 2Ni, < A^r < 4A^l, ^ffqm < for all 6, but gives 
a non-monotonic magnetoconductance, with a minimum 
at b^ = {Nn~2Ni^)/{2Ni^). 

(c) For A^R > 4A'l, ^,9qm > at small 6, but becomes 
negative for 6^ > (A^r, - 4A^l)/(4A^l), and then goes to 
zero for large b. As in (b), the curve is non-monotonic 
with minima at b"^ = (A'r - 2A^l)/(2A^l)- 
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Figure 3: Magnetoconductance curves for the set-up of 
Fig.[T]3. Left panel: fceT = O.IE^t, and iVa/iVL = n + 0.2, n = 
0, 1,2, ...7 (from bottom to top). Right panel: Ny^INi, = 0.2 
(dashed red) and 7.2 (dot-dashed blue), for kTiT/E-v = 0.1, 1, 
2, 4 and 8 (dashed: from bottom to top; dot-dashed: from top 
to bottom). For both panels, the vertical axis gives Sg^m in 
units of the conductance quantum 2e^//i with channel num- 
bers chosen such that NgNl/{N-L + N-r.)^ = 5 in all instances. 
Note the crossover from monotonic to non-monotonic behav- 
ior of the magnetoconductance as T increases, for Nr/Nj^ = 
0.2 (dashed red curves). 



These different regimes persist at finite temperature as is 
illustrated in Fig. [3l however, the boundary between (a) 
and (b), as well as the positions of the minima of the mag- 
netoconductance curves arc T-dependent. The conduc- 
tance exhibits a metal-like (insulating-like) behavior in 
the form of a decrease (increase) of the conductance with 
T, depending on the sign of (Ar — 4(l-ffe)'^AL). This sign 
can easily be changed, whenever one has control over the 
lead widths or the magnetic flux. Remarkably, a mono- 
tonic magnetoconductance may become non-monotonic 
upon increase of the temperature (dashed red curves in 
Fig. [3]). 

Set-up with an S island. In the second of the two 
possible two-terminal set-ups, Fig.[TJ;), the quantum dot 
is connected to a superconducting island, whose chemical 
potential is floating, and adapts itself to a value guaran- 
teeing current conservation, /l = —Ir- Using the ex- 
pression in Ref. (3] for the two-terminal conductance in 
terms of the transmission probabilities, , we obtain 
9 = + Sg'^UT, b) where gjf = NM{N^ + A^r) and 
Km = A^lA^rA^|(^l + A^r)^' f{T^)- This reproduces 
the random matrix theory result [10[ to leading order in 
[A's/(Al + Ar)]^. This quantum correction always in- 
crease the conductance (antilocalization) by a paramet- 
rically large amount (many e^/ft-), with a monotonic mag- 
netoconductance curve. 

Mesoscopic conductance fluctuations and cur- 
rent noise. Ref. [1^ used random matrix theory to 
show that conductance fluctuations remain 0(e^//i) in 
the presence of superconductivity. Our Feynman rules 
reproduce this result. Contributions to var[(;] are the 
product of any two contributions in Fig. [5] connected by 
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encounters. Since path-pairs are not swapped at entrance 
and exit, the connection must involve at least two addi- 
tional encounters with four additional path-pairs, and 
the resulting contribution behaves as N^'^ times the av- 
erage conductance squared. This is at most 0[iV:^], thus 
the quantum corrections to the average conductance are 
parametrically larger than the conductance fluctuations, 
and are therefore easily observable. 

The S contact also leads to e-h contributions to the 
current-noise [ISl . The Feynman rules show that they are 
0[7Vt(^s/^t)"] forn > 1 and are thus smaller than the 
0[A^t] e-e contributions which give the noise in the ab- 
sence of an S lead. Therefore, to leading order in iVg/iYr, 
the parametric magnitude of the zero-frequency current- 
noise is unaltered by the presence of the S lead. 

Effect of SOI. Spin-orbit interaction (SOI) can be 
treated as rotating the spin along otherwise unchanged 
classical trajectories, multiplying Eq. (pi) by Tr[d^ic?^2]j 
where d-yi is the SU(2)-phase of path 71 [25[. For eeO, ee2I 
and hel, this gives a factor of two for spin-degeneracy, 
because d^i = d^2- However for ee2II and he2III it gives 
Tr[c?f], and for hc2I and hc2II it gives Tr[c?ffi2]: where 
di, c?2 are statistically independent random SU(2) phases. 
When the SOI time is shorter than td , one averages uni- 
formly over the SU(2) phases [2^, which multiplies ec2II 
and he2III by -1/2, and he2I and he2II by 1/4. Taking 
T = and neglecting Sg^ii for simplicity, we find that 
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TVl [(1 - 2/(3)^Nn + 4(1 - 2//3)N^] Nj 

(A^L + iVR)3 

(1 - 2/f3)^NLNRNi/{NL + Nnf, 
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for the three standard symmetry classes, /? = 1 (TR sym- 
metry without SOI); 2 (no TR symmetry) and 4 (TR 
symmetry with SOI). Note the presence of the same sym- 
metry prefactor (1 — 2//3) as for weak localization with- 
out superconductivity. Thus with SOI (/3 = 4), both 
Sgqm and (5(7q^ always enhance conductance. Therefore, 
SOI must be absent for a sign change of Sg^m with lead 
width. Turning on SOI (going from /3 = 1 to /3 = 4) 
never changes the sign of (5(7q,',j but changes the sign of 
5gq^ for Nn < 47Vl. 

Concluding remarks. The derivation outlined here 
is for ballistic quantum dots, however the Feynman rules 
that we analyze apply to any system well-modelled by 
random matrix theory. Thus our results are equally ap- 
plicable to disordered dots. We also expect qualitatively 
similar behaviors in diffusive metals coupled to supercon- 
ductors at intermediate temperatures, fceT ^ E^- Work 
in this regime is in progress. 

Upon completion of this work, we noted Ref. 
which uses a somewhat similar methodology as ours in 
closed Andreev billiards. RW thanks L. Saminadayar and 
C. Bauerle for stimulating discussions, and access to their 
data 0. PJ thanks the Physics Department of the Uni- 
versities of Geneva and Basel as well as the Aspen Center 
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